
/* 
robustness to additional outcomes
*/

set more off

estimates clear
clear

use ".\Dropbox\Vaccination\covid.dta" // load the data

drop if statefips==48|statefips==15 // drop states without the vaccination data

zscore(rep*2016 vacc* socialdist_avg* mask hesitant) // normalize the variables

label variable z_rep_pres2016 "Trump vote share, 2016"
label variable z_rep_house2016 "Republican vote share, house 2016"

local county "lat lon temp rain income native_share white_share col_share elder_share male_share  mfg_share_diff popden" // socio-economic controls


// OLS, Trump vote share

 reg z_hesitant z_rep_pres2016 i.statefips, cluster(statefips) 
estimates store hesitancy_geo

 reg z_hesitant z_rep_pres2016  `county' i.statefips, cluster(statefips) 
estimates store hesitancy_all

 reg z_mask z_rep_pres2016 i.statefips, cluster(statefips) 
estimates store mask_geo

 reg z_mask z_rep_pres2016 `county' i.statefips, cluster(statefips) 
estimates store mask_all

 reg z_socialdist_avg z_rep_pres2016 i.statefips, cluster(statefips) 
estimates store socialdist_geo

 reg z_socialdist_avg z_rep_pres2016 `county' i.statefips, cluster(statefips) 
estimates store socialdist_all


estout hesit* mask* social* using ".\Dropbox\Vaccination\Draft\tab_alternative_ols.tex", replace style(tex) cells(b(fmt(3) star) se(par fmt(3))) stats(r2 N, fmt(2 0) labels("R-squared" "Observations")) keep(z_rep_pres2016) label mlabels(none) collabels(none) starlevels(* 0.10 ** 0.05 *** 0.01)

